'''
Created on Jan 29, 2012

@author: nzhao
'''
import math
import numpy as np

Angstrom = 1.0 * 10**-10;
Mu0 = 4 * math.pi * 10**-7;
Hbar = 1.054571596 * 10**-34;

NV_ZFS = 2.87 # GHz
NV_N14_ZFS = -5.04 # MHz;

# gyromagnetic ratios in rad / s / T
GAMMA_E = -1.760859794 * 10**11 # free electron 
Gamma_C13 = 6.726149871 * 10**7;# Carbon-13
Gamma_N14 = 1.9338 * 10**7 # Nitrogen-14
Gamma_N15 = -2.7116 * 10**7 # Nitrogen-15
Gamma_H = 2.67522212 * 10**8 #Proton

X=np.array([[0.,1,],[1.,0.]])
Y=np.array([[0.,-1.0j],[1.j,0.]])
Z=np.array([[1.,0,],[0.,1.]])

SPIN1X = np.array([ [0., 1., 0.], [1., 0., 1.], [0., 1., 0.] ]) / math.sqrt(2.)
SPIN1Y = 1.j * np.array([ [0., -1., 0.], [1., 0., -1.], [0., 1., 0.] ]) / math.sqrt(2.) 
SPIN1Z = np.array([ [1., 0., 0.], [0., 0., 0.], [0., 0., -1.] ])

X_DIR = np.array([1., -1., 0.]) / math.sqrt(2.)
Y_DIR = np.array([1., 1., -2.]) / math.sqrt(6.)
Z_DIR = np.array([1., 1., 1.])  / math.sqrt(3.)
